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Cq ' Abstract 

The process e~^e~ — > Htt can be used at the Next Linear Colhder to measure the 
Higgs-top Yukawa couphng. In this paper, we compute 2 — > 8 processes of the form 
e"'"e~ —>■ bbbbW~^W~ — > bbbbl^Vi qq', accounting for the Higgs-top-antitop signal as 
well as several irreducible backgrounds in the semi-leptonic top-antitop decay channel. 
We restrict ourselves to the case of a light Higgs boson in the range 100 GeV ~ Mh ~ 
140 GeV. We use helicity amplitude techniques to compute exactly such processes at 
tree level in the framework of the Standard Model. Total rates and differential spectra 
of phenomenological interest are given and discussed. 
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At the Next Linear Collider (NLC), running with a centre-of-mass (CM) energy of ^/s = 
500 GeV [m, the Higgs boson of the Standard Model (SM) can be produced in association 
with top-antitop pairs [0, through the process e~^e~ —>■ Hti, which proceeds via the diagrams 
displayed in Fig. |l[ That is, the scalar particle can be radiated either from the top quark 
pair or from a virtual Z boson. In the latter case, it is the neutral gauge vector to eventually 
produce the heavy quark pair. Clearly, given the actual value of the top mass, mj ^ 175 
GeV, between the two sets of graphs, it is the first one which dominates. On the one hand, 
the Z* — * ti decay occurs far off the mass-shell of the Z boson. On the other hand, the large 
Yukawa coupling exceeds the strength of the HZZ vertex. Indeed, it is the possibility of 
measuring such Yukawa interaction that renders associated production of Higgs bosons and 
top (anti)quarks phenomenologically interesting at the NLC 0. 

From the above values of y/s and rrit, it follows that only Higgs scalar s with mass Mh 
up to 140 GeV or so can be produced, because of the kinematical limit imposed by the 
difference y^ — 2mt. For such values of Mh, the dominant Higgs decay mode is if — >• bb, 
this being overtaken by the off-shell decay into two ly^'s, i.e., H -^ W~^*W~*, only for 
Mh ~ 130 — 140 GeV, see Fig. 1 of Ref. [Q. However, these Higgs masses are extremely close 
to the kinematical limit of the Hti intermediate state, so that the production cross section of 
the latter is very small |^. Furthermore, notice that in order to reconstruct the Higgs mass 
one would require a fully hadronic decay of the W^*W~* pairs produced in the Higgs decay, 
this leading to a signature with at least eight jets in the final state. In fact, at least one top 
quark would be required to decay into jets, in order to exploit the reconstruction of its mass 
to reduce various QCD backgrounds. In other terms, the search for H -^ W~^*W~* decays 
from e'^e" -^ Hti would be of difficult experimental use, considering the reduced number 
of events, the rather chaotic topology and the problem that the latter generates, because of 
the combinatorics, while attempting to disentangle the H and t resonances. In the end, one 
would be much better off to rely on the two-body mode H —>■ bb over the entire Mh range 
allowed by Higgs-top-antitop intermediate states at i/i = 500 GeV. 

As for ti decays, one would most likely exploit the semi-leptonic channel, i.e., ti — > 
bbW^W~ — > bbi^Ui qq', where i and z/ represent a lepton at high transverse momentum (to 
be used for triggering purposes) and its companion neutrino and qq' refers to the two possible 
combinations of light quark pairs and their charge conjugated channels (neglecting Cabibbo- 
Kobayashi-Maskawa mixing effects). This is the decay signature we will concentrate on. As 
a matter of fact, such a choice is not restrictive, in the sense that the latter is to date the 
experimentally preferred channel in searching for ti ^ bbW^W~ events P]. 

If one does assume such Higgs and top decay modes, then signal events can be searched 
for in data samples made up by four b quark jets, two light quark jets, a lepton and a 
neutrino. In other terms, a '46 + 2 jets + £^ + -Emiss' signal, assuming the four heavy quark 
jets to be recognised as such thanks to the /x-vertex devices of the NLC, with £ = e, /i, r| 
and where the missing energy, -Emiss, originates from the neutrino escaping detection. 

Though the calculation of the on-shell production e^e~ —>■ Hti has been tackled long ago 
[^ , that of the complete 2-^8 body reaction, without any factorisation of production and 
decay processes, has never been attempted before. Not surprisingly so, as even in presence of 
only five diagrams, both the large number of particles in the final state and the complicate 

^We include r's to enhance the signal rate, assuming that they are distinguishable from quark jets. 



resonant structure of the latter impose non-trivial problems to the matrix element (ME) 
calculation and to its integration over the phase space, respectively. Things become even 
more involved if one starts including (irreducible) backgrounds in the calculation, as needed 
in order to realistically simulate phenomenological studies. For example, if one restricts 
oneself to all those channels that proceed through an intermediate HbbW^W~ stage, then 
the full gauge invariant set (including Higgs bosons produced via other graphs than those in 
Fig. |l|) counts 350 tree-level diagrams. 




Figure 1: Relevant Feynman diagrams contributing at lowest order to the process 6^62 — ^ 
HstJ:^. An internal wavy line represents a 7 (graphs 1, 3) or a Z (graphs 2, 4, 5). 

One of the main irreducible backgrounds to the Higgs-top-antitop signal at the NLC 
is the scattering e^e~ -^ Zti [^, if one considers that the two processes have comparable 
production cross sections P, H and that the Z boson decays into hh pairs some 15% of 
the times. Even though the difference between Mh and Mz is always larger than 10 GeV 
(assuming a late 100 GeV bound on the former from the all of LEP2 data 0) and the 
width of the Higgs boson is very narrow (about ten MeV at the most for masses up to 
140 GeV, see Fig. 2 of Ref. [§), one should recall both the large value of that of the Z 
boson, Tz ~ 2.5 GeV, the finite efficiency of the detectors in reconstructing jet energies and 
directions (to say the least, yielding a resolution of some 5 GeV in invariant mass) and the 
mis-assignment problems arising when pairing the four h jets in the final state in the attempt 
to recognise resonances in the hh decay channel. Thus, it is inevitable to conclude that Zti 
events will represent a serious noise. On-shell Z-top-antitop production proceeds at tree- 
level through the nine graphs of Fig. 0. If one however considers, on the same footing as was 
done for Higgs production, all the gauge invariant set of amplitudes producing ZbbW^W~ 
intermediate states, followed by Z ^ 66, then the number of graphs involved is 546. (Notice 
that several of the production channels described by the latter do involve Higgs bosons, some 
of which decay into hh pairs.) 

In addition, one should also consider e^e" -^ ghhW^W~ intermediate states, where g 
represents a gluon eventually yielding hh pairs. Although none of hh invariant masses has 
in this case the tendency of being produced around Mh (in particular, the one induced by 
the g splitting logarithmically increases at very low mass values, because of the infrared 
singularity of QCD, only regulated by the h mass, m^), such mechanisms proceed through 
strong interactions, so that their production rates could well be comparable to those of the 
signal^]. In fact, because of the mis-pairings of h quarks, large tails in the hh invariant mass 
distributions could arise, despite of the softness and collinearity of two of the heavy quarks. 
The dominant background contribution from these mechanisms would come from e+e~ — > gti 
events |T^, with the gluon radiated before the (anti)top decays take place. There are four 

■^For opposite reasons, one can avoid studying e^e~ -^ jbbW~^W~ reactions, with the photon sphtting 
into bb pairs. 



tree-level diagrams associated with this 2-^3 process, see Fig. ^ The total number of those 
yielding gbbW^W~ states is instead 152. 




Figure 2: Relevant Feynman diagrams contributing at lowest order to the process efe2 — > 
Z^t/ii^. An internal wavy line represents a 7 (graphs 1, 3, 5, 7) or a Z' (graphs 2, 4, 6, 8, 9). 







Figure 3: Relevant Feynman diagrams contributing at lowest order to the process ei'e2 
Q'itii^. An internal wavy line represents a 7 (graphs 1, 3) or a Z (graphs 2, 4). 



It is the purpose of this letter to compute all such processes and compare the signal rates 
and distributions to those obtained from the various backgrounds that we have described, in 
order to assess the chances of genuinely exploiting the Higgs-top-antitop production process 
in measuring the Higgs-top Yukawa coupling. In this respect, the reader should notice one 
subtlety. In fact, the mentioned coupling not only appears in the e~^e~Hti — >■ HbbW~^W~ 
'signal', but also in several 'background' mechanisms, such as in e+e" -^ W'^*W'^ production 
with one of bosons off-shell, followed by W^* -^ Htb + Hib -^ HbbW^ and e+e" — > t*bW~ 
production of an off-shell t quark, eventually yielding t* -^ Ht —>■ HbW^ (plus the charged 
conjugate case). These can be regarded as 'single top' processes, as opposed to the 'double 
top' one, i.e., e~^e~ — > Hti, themselves being proportional to the Higgs-top Yukawa coupling. 
More correctly then, these two subprocesses should be considered as additional contributions 
to the, say, 'Yukawa' signal, further recalling that they carry one resonant top decay (we 
are selecting the semi-leptonic channel, thus implicitly assuming that no more than one top 
mass can in principle be reconstructed). 

To compute all signal^ and background graphs we have resorted to helicity amplitudes 
methods. In particular, we have made use of the HELAS subroutines [|1^, based on the 
formalism of Ref. |13|. All the FORTRAN codes produced this way have been tested for gauge 
invariance satisfactorily, so to give us confidence in our numerical results. Furthermore, 

*Note that we calculate the Higgs-top-antitop signal at the leading-order (LO), though we are aware that 
several higher order corrections (mainly to the on-shell production) are known to date [g, |l^ . We do this 
for consistency, as all the 2^8 background processes are evaluated here at tree level. 



the 2^8 'dominant' (as we shall see below) signal and background processes of the form 
e'^e^ -^ XbbW~^W~ -^ bbbbi^i'iqq' , with X = H, Z and g, have also been implemented 
by using the spinor techniques described in Refs. [H, 15]. Wherever the two approaches 
overlapped, we have seen perfect agreement between the outputs of the two sets of codes. 

Numerical results have been produced after integration of the Feynman amplitudes 
squared over eight-body phase spaces. In order to account accurately for all their com- 
ponents, we have split the MEs of the form e"'"e~ —>■ XbbW~^W~ -^ bbbbi^uigq' in resonant 
sub-terms and integrated each of these separately. Only in the end the various integrals were 
summed up, in order to recover gauge-invariance |jl6|. The algorithms used to perform the 



multi-dimensional integrations were VEGAS [^ and, for comparison, RAMBO [T^ 



To describe the vector and axial couplings of the gauge bosons to the fermions, we have 
used sin^ 9w = 0.2320. The strong coupling constant a^ entering the QCD processes (i.e., 
X = g) has been evaluated at two loops, with Nf = 4 and A^jg = 230 MeV, at a scale 
equal to the collider CM energy, y/s = -Eecm = 500 GeV. The electromagnetic coupling was 
Oem = 1/128. For masses and widths, we have used: m£ = rrii,^ = rriu = rrid = iris = rric = 0, 
rrib = 4.25 GeV,mt = 175 GeV (as default), Mz = 91.19 GeV, F^ = 2.50 GeV, Mw = 
80.23 GeV and F^y = 2.08 GeV. As for the top width F^, we have used the LO value of 1.5 
GeV. Only in one circumstance, in order to study the sensitivity of the signal processes to 
the Higgs-top Yukawa coupling, we have changed rrit by ±5 GeV. The widths corresponding 
to these two new values are 1.3 and 1.6, for the lower and higher rrit figure, respectively. 

Concerning the Higgs boson, we have spanned its mass Mh over the range 100 to 140 
GeV. As for its width, Th, we have computed it by means of the same program described 
in Ref. [Q, which uses a running b mass in evaluating the H ^ bb decay fraction. Thus, for 
consistency, we have evolved here the value of rrih entering the Hbb Yukawa coupling of the 
H ^ bb decay current in the same way as then. 

Finally, notice that starting from our 2^8 MEs for e~^e~ —>■ Xti —>■ bbbbW~^W^ -^ 
bbbbi^Uiqq', in all cases X = H, Z and g, we are able to reproduce (apart from minor spin 
correlations) the cross sections that one obtains from the 2 — *■ 3 ones for e~^e~ — ^ Xti, times 
the relevant branching ratios (BRs), by adopting a Narrow Width Approximation (NWA) for 
the various resonances R involved (i.e., R = H,t, W^ and Z), by rewriting the corresponding 
(denominator of the) propagators as (for F = F^ the standard expression is recovered): 

1 / F y/2 

p2 — m| + im/jF VF/j/ 

with F — i> 0, this way mimicking a delta distribution, i.e., (5(p^ —mj^). (In the case X = g we 
had to supplemented the 2^3 ME for e~^e~ — ^ gti with the splitting function for g ^ bb.) 

In the following, total and differential rates are those at parton level, as we identify jets 
with the partons from which they originate. Gaussian smearing effects are simulated. No 
efficiency to tag four b quarks is included. 

We start our analysis of the results with a disclaimer: we have not included Initial State 
Radiation (ISR) [|19| in our calculations. We have done so mainly for technical reasons. 
Simply because we are already dealing with complicated processes requiring delicate inte- 
grations, over nineteen dimensions and with a laborious rearrangement of the phase space, 
to account for the multi-resonant behaviour of hundred of diagrams, that even adding the 



ISR in the simplest wayQ would prove rather costly in terms of efficiency of the computa- 
tion. In addition, we would expect ISR to affect rather similarly the various processes of 
the form e"'"e~ — » XbbW~^W~ — > bbbbi^Uigq'. As we are basically interested in relative rates 
among the latter, we are confident that the basic features of our results are indifferent to 
the presence or not of photons radiated by the incoming electron-positron beamsQ. 
Fig. H presents the production cross sections for the following (sub)processes: 

1. ffistly, the 2^3 on-shell ones, 

e+e' -^ mi, (2) 

e+e" -^ Zti, (3) 

e+e- ^ gti, (4) 

as obtained from the diagrams in Fig. |^-^ multiplied by the BRs and the g —^ bb 
splitting function; 

2. secondly, the 2^8 ones which proceed via those above, 

e+e' -^ Hti -^ HbbW+W- -^ bbbU+vm', (5) 

e+e" -^ Ztt -^ ZbbW+W- -^ bbbbtvm\ (6) 

e+e- -^ gti -^ gbbW+W' -^ bbbU+uggq', (7) 

as obtained from the diagrams in Fig. |]-|^ supplemented with the decay currents; 

3. thirdly, the 2^8 ones including also all other diagrams, 

e+e- -^ bbW+W- -> bbbU+vm', (8) 

e+e- -^ ZbbW+W- -^ bbbU+ueqq', (9) 

e+e- -^ gbbW+W- -^ bbbW+vm' ■ (10) 

In moving from cases 1. to 2., one can appreciate the onset of spin and width effects, see 
top of Fig. §, whereas in comparing 2. and 3. one can disentangle those due to the diagrams 
not proceeding via X-top-antitop pairs, see bottom of Fig. ^. It turns out that spin and 
width effects are sizable only for the Higgs-top-antitop and Z-top-antitop processes, not for 
the gluon-top-antitop ones. They are of the order of +&% in Zti diagrams, whereas in the 
case of Hti they vary between +2.5% at Mh = 100 GeV and -15% at Mh = 140 GeV. For 
gti diagrams they amount to less than 1% (hence the overlapping of the two dotted curves 
in the top frame of Fig. ^). As for effects due to non-X-top-antitop graphs, things go the 
other way around. The gti rates are hugely increased, by as much as a factor of two, whereas 
the Zti and Hti ones never get larger than 2.3% and 4.3%, respectively. The growth of the 
QCD rates is mainly due to the large amount of gluon radiation (here, eventually yielding 
bb pairs) produced in the top quark decays [|10|- Note, however, that the latter can easily 



be controlled by imposing that none of the invariant masses of five particle systems with 



^For example, via the so-called Electron Structure Function (ESF) approach [ |l9| . 

^We also neglect beamsstrahlung and Linac energy spread, by assuming a narrow beam design ||l£ 



three heavy and two hght quarks (and/or two leptons, if the vi niomentum is reconstructed) 
reproduces m^. 



Therefore, by studying the production rates of all reactions (0)-(lO), one may remark on 
two key aspects. On the one hand, the bulk of the cross sections of processes (|8|)- (p!oD comes 
from the X-top-antitop channels (0)-(|3)- On the other hand, the QCD process (0) is the 
dominant one, for any value of Mh. (That for M^ ~ 125 GeV or so the Ztt production rates 
started exceeding the Hti ones was rather trivial to derive |@, H].) Whereas the first result 
was clearly expected, the second one came as somewhat of a surprise. As a consequence, 
in the reminder of our analysis, we will mainly concentrate on the X-top-antitop diagrams 
and study some of their differential spectra that can help disentangling the Higgs diagrams 
from the Z and, especially, the gluon ones. We will do so for the choice Mh = 130 GeV, as 
representative of the case in which both Z-top-antitop and gluon-top-antitop backgrounds 
overwhelm the Higgs-top-antitop signal (see Fig. ^). 

As we have already stressed that one of the hb pairs in the final state would naturally 
resonate at Mh, at Mz or logarithmically increase at low mass, for processes (^), (H) and (^, 
respectively, we start investigating the di-jet mass spectra that can be reconstructed from 
the four h quarks in the '46 + 2 jets + £^ + -Emiss' signature. Since we do not assume any jet- 
charge determination of the (/i-vertex tagged) h jets and consider negligible the mis-tagging 
of light-quark jets as heavy ones, six such combinations can be built up. We distinguish 
among these by ordering the four h jets in energy (i.e., Ei > E2 > E3 > E4), in such a 
way that the 26 invariant mass rriij refers to the ij pair (with i < j = 2,3, 4) in which 
the i-th and j-th most energetic particles enter. Having done so, one should expect to see 
the typical resonant/logarithmic behaviours described above now 'diluted' in the various ij 
combinations. This is evident from Fig. ^. There, one can appreciate the resonant shapes 
around Mh and Mz in all ij cases (for ij = 12, the Z peak is just a tiny kink on top of a 
Jacobian shape). As for the 'divergence' in the g ^ bb splitting of the QCD process, this can 
easily be spotted in the case ij = 34. In the end, the 26 mass spectra look rather promising 
as a mean of reducing both backgrounds (|D-(|3)- By requiring, e.g., mu > 50 GeV, one 
would vigorously reduce the latter; similarly, by imposing, e.g., |mi4 — Mz\ > 15 GeV one 
would reject the former considerably. 

Another way of looking at the same phenomenology in processes (|^)-(0) is by studying 
the energy spectra of the four 6 quarks. In fact, the larger value of Mh, as compared to 
Mz, should boost the 6 quarks generated by the Higgs boson towards energies higher than 
those achieved in the Z decays. Conversely, the energy of the 6 quarks emerging from 
the two remaining unstable particles, top and antitop quarks, should be softer in the first 
case. Following similar arguments, one should expect the hardest (softest) 6 (anti)quark from 
gluon events to actually be the hardest (softest) of all cases (|^)-(0), once again, because of 
the infrared QCD splitting of a soft gluon. Recalling that the two most energetic 6's seldom 
come from a. H, Z or g splitting in XbbW~^W~ intermediate states (see top-left curves in 
Fig. ^), the above kinematic features are clearly recognisable in Fig. ^. Therefore, the energy 
spectra too are rather useful in disentangling Higgs events. If one imposes, e.g., Ei < 100 
GeV and E^ > 50 GeV, both Z and gluon events can be strongly depleted, at a rather low 
cost for the signal. 

We conclude the numerical analysis by studying the sensitivity of the signal to the Higgs- 
top Yukawa coupling, by varying the top mass by 5 GeV above and below its default value. 
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0-tot (fb) 


mt (GeV) 


HhhW+W' 


ZhhW+W- 


ghhW+W- 


170 


0.061 


0.087 


0.79 


175 


0.033 


0.067 


0.77 


180 


0.012 


0.050 


0.74 



Table 1: Cross sections for processes (see the text) (|])-(0), for three discrete values of the 
top mass. The total CM energy is ^/s = Ecm = 500 GeV. The Higgs boson mass is Mh = 130 
GeV. No cuts have been implemented. (Numerical errors from the Monte Carlo integration 
do not affect the significant digits shown.) 

However, as to modify rrit (and, consequently, Fj) has also incidence on the top propagators, 
and since these enter many of the diagrams associated with processes (^- (p^Of ), we present 
the rates for the latter, that is, for the full sets of diagrams in each case. This is done in 
Tab. |l|. The value we have chosen for the Higgs mass, i.e., Mh = 130 GeV, is a critical one 
for process (^, the main source of events (|^). In the sense that the sum 2mt + Mh is very 
close to y/s, so that the corresponding rates in Tab. |l] (see second column) are the result of 
the interplay between the rise of the cross section with mf and its fall because of the phase 
space suppression (width effects are less relevant). Indeed, between the two tendencies is 
the latter to dominate. In fact, the production cross sections of all three processes (H)- 
( p!OD decrease with increasing top mass. Even in presence of such delicate interplay, the 
sensitivity of Higgs-top-antitop events to the actual value of the top mass is rather strong, 
as the corresponding cross section changes by a factor of 5 between rrit = 170 and 180 GeV. 
Backgrounds variations are always smaller. However, both processes (§)-([l3) are larger than 
(H). Once again, it has to be stressed that background rates ought to be reduced severely if 
one wants to perform dedicated studies of the Higgs-top Yukawa coupling. 

In summary, in our opinion, the study of the Higgs-top Yukawa coupling at future 
electron-positron colliders, such as the NLC running with a CM energy of 500 GeV, can 
in principle be pursued by means of the Higgs-strahlung process e~^e~ — > Hti. In fact, the 
irreducible backgrounds affecting the latter can be brought under control in the semi-leptonic 
top-antitop decay channel ti -^ bbW~^W~ -^ bbC.^i'^qq' , further assuming H -^ bb, as natural 
for Higgs masses up to 140 GeV or so. 

However, this requires to somehow recognise the b jets in the final state with high effi- 
ciency, as the observable rates of the signal are below the femtobarn level. The knowledge of 
the momenta of the heavy quark jets entering the signature '46-1-2 jets -t-£^-|-i5miss' is crucial 
in order to reduce the overwhelming QCD background, mainly proceeding via e~^e~ -^ gti 
events, if one aims to disentangle such Higgs events at all. The competing electroweak back- 
ground, mainly proceeding through e+e^ -^ Zti intermediate states, can be dealt with if 
the mass resolution of di-jet pairs of b quarks is around 10 GeV or less. Other irreducible 
background channels, induced by e'^e~ -/^ Xti -^ Xbb W^W~ intermediate states, with 
X = H, Z or g, are significantly smaller than those proceeding via X-top-antitop graphs, 
with the only exception of QCD graphs involving one radiative (anti)top decay. 

In the end then, although a careful simulation of possible tagging strategies should even- 
tually be performed, we believe that, if the Higgs mass turns out to be in the intermediate 
range, the NLC constitutes an ideal laboratory for the kind of studies sketched here. We 



base our conviction on the fact the we have performed a new and rather complete calculation 
of signal and backgrounds involving up to ten external particles. 
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Figure 4: Cross sections for the following processes (see the text): (||) (long-dashed line, 
above), (H) (dot-dashed line, above) and @) (dotted line, above); (||) (solid line, above and 
below), (H) (dashed hne, above and below) and (|^ (fine-dotted line, above and below); (|D 
(long-dashed line, below), (|^) (dot-dashed line, below) and (|T0D (dotted line, below). The 
total CM energy is ^/s = Ecm = 500 GeV. No cuts have been implemented. 
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Figure 5: Differential distributions in the invariant mass of all possible combinations ij, with 
i < j = 2, ...4, of the energy-ordered h jets (i.e., such that Ei > E2 > E^ > E4) for the 
following processes (see the text): (^ (solid line), (H) (dashed line) and (^ (dotted hue). 
The total CM energy is ^/s = Ecm = 500 GeV. No cuts have been implemented. Note that 
the rates of reaction (|^ have been divided by three for readability. 
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Figure 6: Differential distributions in energy of tlie energy-ordered b jets (i.e., sucli tliat 
El > E2 > E^ > E4) for tlie following processes (see the text): (||) (solid line), (||) (dashed 
line) and (^ (dotted line). The total CM energy is y/s = E^m = 500 GeV. No cuts have been 
implemented. Note that the rates of reaction (|^) have been divided by three for readability. 
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